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There were two major developments in 1994 in our collaborative GPS experiment in 
the Tien Shan of the Former Soviet Union (FSU). Both were motivated by our expecta- 
tion that we will ultimately obtain better science at lower cost if we involve our col- 
leagues in the FSU more deeply in 1) the collection and 2) the analysis of data. As an 
experimental test of the concept of having our local collaborators carry out the field work 
semi-autonomously, we sent 6 MIT receivers to the Tien Shan for a period of 3 months. 
To enable our collaborators to have the capability for data analysis, we provided comput- 
ers for two data analysis centers and organized a two-week training session, run by 
Professor Thomas Herring of MIT. 

This report emphasizes the rationale for deeper involvement of FSU scientists, 
describes the training sessions, discusses the data collection, and presents the results. We 
also discuss future plans. More detailed discussion of background, general scientific 
objectives, discussions with collaborators, and results for the campaigns in 1992 and 
1993 have been given in previous reports. This project involves the amalgamation of 
two, separately funded projects, one by Pis Hamburger and Reilinger and the other by Pis 
Molnar and Hager. Each of these two projects is funded jointly by NASA and NSF. In 
addition, the work is being conducted under the auspices of the US-USSR Agreement on 
Cooperation in the Field of Environmental Protection, with support from the United 
States Geological Survey. 
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GENERAL SCIENTIFIC OBJECTIVES AND BACKGROUND 

As the world's most prominent intracontinental mountain belt, the Tien Shan provides 
an outstanding natural laboratory to address two general problems in geodynamics: (1) 
the processes by which distributed intracontinental shortening occurs, and (2) the mecha- 
nisms that allow strain to be absorbed adjacent to an intracontinental strike-slip fault 
along which slip varies markedly. Moreover, with a high level of seismicity and with 
rapid deformation, the Tien Shan offer the potential for witnessing post-seismic deforma- 
tion and the opportunity to study localized deformation within the regional framework of 
intracontinental tectonics. 

Intracontinental mountain building. The present range owes its height and high level 
of seismic and tectonic activity to roughly north-south crustal shortening. Hence, it 
serves as a prominent consequence of India's collision with and subsequent penetration 
into the Eurasian Plate (Figure 1). Yet, this area lies 1000-2000 km from the northern 
edge of the Indian plate, and hence deformation should be largely unaffected by peculiar- 
ities of the ancient plate margin. Thus, the Tien Shan provides a laboratory to study 
intracontinental shortening without the influence of asymmetries of kinematics and 
structure established at the time of the collision, as are clear, for example, in the 
Himalayas and the Alps. 

Approximately two thirds of our network spans a segment of the Tien Shan where the 
structure is relatively simple. Ranges and intervening basins trend nearly east-west. 
Blocks of crust, with east- west dimensions of approximately 100 km or more, constitute 
ranges that are thrust over other blocks to form basins. North-south dimensions of such 
blocks of only 20-30 km are comparable to the depths to the more ductile lower crust. 
We anticipate that a comprehensive study of the velocity field in this region will allow us 
to study the interconnections of the block movements and deformation at depth and to 
understand how the penetration of India into the rest of Eurasia is communicated to the 
areas farther north. 

Strike-slip geodynamics. The Tien Shan offers a good example of a second feature 
peculiar to intracontinental deformation. A major strike-slip fault, the Talaso-Ferghana 
fault, truncates the western part of the region (Figure 1). The east-west trending blocks 
and the faults between them seem to be cut off by the Talaso-Ferghana fault. Moreover, 
the Talaso-Ferghana fault ends in the northwestern part of the Tien Shan. In oceanic 
regions, transform faults, along which slip is constant, abruptly truncate at spreading 
centers and subduction zones. In continental regions, however, strike-slip faults do not 
end abruptly, but seem to "die out" by distributed, permanent deformation on one or both 
sides of the fault. The Talaso-Ferghana fault is perhaps the world’s premier example of 
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such an intracontinental strike-slip fault. Slip seems to be rapid, possibly 10 mm/yr. Yet, 
the fault passes from being a major structure to being absent in a distance of only 100 km. 
Bounding the fault on both sides are high ranges suggestive of rapid crustal shortening, 
but also possibly reflecting rapid rotations about vertical axes. The geometry of defor- 
mation permits us to determine the kinematics of deformation using GPS. This kinemat- 
ics should facilitate an understanding of the importance of the relative strengths of the 
brittle layer at the surface and the region below it that flows ductilely. 

We began field work in July, 1992, with a small-scale pilot experiment in the north- 
ern Tien Shan, in the republics of Kazakhstan and Kyrgyzstan. In June - July, 1993, we 
carried out a 7-week field experiment which constituted a major expansion of this pilot 
network. The field program included reconnaissance and installation of benchmarks for 
regional GPS sites in Kazakhstan and Kyrgyzstan, a training program for collaborating 
scientists, and an 18-day observation campaign covering the newly installed 86-station 
geodetic network. The network is closely tied to a regional GPS network established by 
the German GeoForschungZentrum geodetic group (GFZ), under the direction of Christof 
Reigber, and Kazakh, Kyrgyz, Russian, and Uzbek scientists. 

As reported previously (Souter et al, 1993), we determined apparent site displace- 
ments for 12 sites by comparing the 1992 and 1993 surveys. To summarize these results, 
we showed that we could determine relative positions between sites hundreds of kilome- 
ters apart quickly, efficiently, and precisely, with formal errors for daily solutions of 
approximately i 5 mm. Observing for many sessions decreases the uncertainties. But 
this decrease is not as rapid as it would be if the error spectrum were white. In particular, 
the largest apparent displacement between the initial 1992 survey and the repeat survey in 
1993 was between two sites located on relatively stable portions of the Kazakh platform: 
Kurdai (KURD) and Azok (AZOK). Taken at face value, the difference suggests that 
KURD moved east ~ 15 mm with respect to AZOK in the interval between the measure- 
ments. The difference between relative positions of 15 mm is large, but it is smaller than 
the total range observed for a shorter line between MEDE and AZOK that was observed 
many times in 1992. With only 3 (consecutive) days of observation at KURD (and most 
other sites) each year, we could not reject the hypothesis that the apparent displacement is 
an unfortunate artifact. In fact, this is our preferred explanation - we expect there to be at 
most very slow (< 1 mm/a) relative movement between these sites, for both sites seem to 
lie on the relatively stable Kazakh platform. 

Given our desire to measure and interpret tectonic displacements, not long-term 
errors, it became clear that it would be preferable to survey our network more frequently 
than once every two years, as originally proposed. But large yearly field campaigns with 
many US participants would be prohibitively expensive. From the large experiment that 
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we ran in 1993, it was clear that there were many local scientists who could carry out 
field measurements as competently as we Pis. Perhaps a better use of resources would be 
to provide GPS receivers for our collaborators to use for longer periods of less intense 
campaigns. In addition, our FSU associates desired to become able to independendy ana- 
lyze GPS data collected in collaboration with us and with others. We therefore planned a 
pilot experiment in 1994, with no additional funding, to determine if this approach would 
work. 

Field Measurements 

In May, 1994, UNAVCO shipped 6 Trimble 4000 SSEs (P/X code) belonging to MIT 
to the Institute of High-Temperature Physics (Russian Academy of Sciences, "IVTRAN," 
under the direction of Yuri Trapeznikov) base near Bishkek, Kyrgyzstan. PI Michael 
Hamburger traveled to IVTRAN to provide training and oversee the planning of the cam- 
paign schedule. After this initial training, the experiment was turned over to 
Trapeznikov's supervision for its duration. (The presence of U.S. scientists during three 
different stages of the field campaign permitted us to provide technical support during the 
experiment and to ensure close communication with American Pis.) 

Data were recorded round-the-clock in 10 60-hour and one 40-hour sessions, carried 
out during a six-week long campaign. Receivers operated at the IVTRAN test site near 
Bishkek throughout the recording period and at Cholponata (CHOL) for all but one of the 
sessions. In addition, several other sites were observed during several consecutive 60- 
hour sessions to provide better assessment of repeatabilities. Over the six-week period, 
we were able to observe at 38 of our 86 sites (Table 1). The GPS receivers operated 
nearly flawlessly for the duration of the experiment. 

The equipment was returned by co-I Bradford Hager at the end of July. It was in 
excellent condition - the only "casualty" being a crashed hard disk in one of the PCs. 
(Hager brought these receivers back to the US as accompanying luggage at a cost of only 
$ 100/container.) 

Data Analysis 

The original data recorded in the field were transferred onto floppy disks. Three 
copies of all field data were made and verified before leaving Bishkek. We left one copy 
with IVTRAN and brought two copies to the US. One copy has been archived at 
UNAVCO. We will also archive a copy with CDDIS. 
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Our data analysis procedure has two main steps. In the first step, we break the rather 
long (60-hour) sessions into a series of 24-hour (or 12-hour) sessions, organized by GPS 
day. We analyze each day of data from all our receivers, along with data from three sta- 
tions in the global tracking network, using the GAMIT software developed at MIT and 
Scripps. We use daily orbits provided by the Scripps orbit facility as the starting model. 
The data are cleaned and cycle slips patched, where possible, automatically. Integer 
ambiguities are resolved, where appropriate, using well-known techniques. The output of 
this stage of the analysis is an estimate of station coordinates and other parameters (e.g., 
orbit perturbations and atmospheric parameters), along with their covariance matrices, for 
each day on which observations were made. 

In the second step of the analysis, we combine the individual daily solutions into a 
solution for coordinates and (where appropriate) velocities using the GLOBK software 
package developed at Smithsonian Astrophysical Observatory and MIT. GLOBK rigor- 
ously combines the adjustments of the site coordinates and orbital parameters and their 
covariance matrices from the individual-day solutions using a Kalman filter to provide 
what we call the "global" solution; we simultaneously estimate local station coordinates, 
global station coordinates, and satellite orbits. 

Our current solution for velocities is shown in Figure 2, where the error ellipses rep- 
resent 95% confidence intervals. The 12 sites in the northern part of the network gener- 
ally have smaller uncertainties because they have been observed 3 times over a two-year 
period, in contrast with the other sites, which have only been observed twice, with one 
year between campaigns. 

For the broad, north-south profile extending across the Tien Shan, the apparent 
deformation rates are sufficiently high (given the presently attainable precision of GPS 
measurements) that they show resolvable displacements within one year of repeat mea- 
surements. Although a second round of repeat measurements is essential to obtain suffi- 
ciently accurate estimates of station position for robust tectonic interpretation (see dis- 
cussion below), the comparison of 1993-94 site positions for the Tien Shan network is 
providing tantalizing initial results that we cannot resist commenting upon: 

(1) The long, north-south baselines extending across the central Tien Shan suggest con- 
sistent north-south shortening, taking place at rates of 15-22 mm/yr. This high rate is 
somewhat surprising, as it accounts for some 30-40% of the total India-Eurasia plate 
convergence rate, in spite of the fact that the network does not span across the 



6 


boundary between the Tarim Basin and the Tien Shan, where much of the present-day 
seismicity is concentrated. 

(2) Most of convergence appears to be accommodated by intramontane shortening within 
the Tien Shan, rather than across the northern foreland of the mountain belt. In par- 
ticular, the largest gradient in the deformation field takes place across the Issyk-Kul 
basin, which separates the Terskei Alatau and the Kungei Alatau ranges in the central 
Tien Shan. 

(3) A systematic divergence of displacement vectors is observed in the southwestern 
portion of the study area, suggesting local block rotations in the southern Tien Shan. 

But these results are preliminary the main focus of this report is on the impact that 
the third year of data has on our estimates of relative site velocities. In Figure 3, we show 
the evolution of our most frequently observed line, between Cholponata (CHOL) and 
Azok (AZOK), where we have had at least 7 days of observations in each of 3 years. The 
daily solutions all cluster about the linear trend. Note, however, the outlier in the north 
component in 1994 - if this line had only been observed on that day, or even on three 
days including that day, the estimate for the northward site velocity would have been sub- 
stantially lower. Also note that the error bars and scatter in the daily solutions in 1994 
are better than for 1992 (primarily due to better orbital control by the global network) and 
only slightly worse than for 1993. Thus, we conclude that we are not paying much of a 
penalty in short-term scatter by observing with a network of 6, rather than 18 instruments, 
a conclusion that has guided our future observing strategy for the Tien Shan network. 


In contrast, the line between Azok and Medeo (MEDE) was observed many times in 
1992, but only 3 times in each of 1993 and 1994 (Figure 4). Based on analysis of the first 
two years of data, we had inferred N-S extension between these two sites — a perhaps sur- 
prising result for an environment dominated by compressional deformation. But the 1911 
M=8.6 Chon Kemin earthquake occurred just to the south of Medeo, and we had provi- 
sionally interpreted this extension as the result of postseismic viscoelastic relaxation in 
the lower crust. The addition of additional data strongly suggests that the inference of 
extension was incorrect. 

As a final example, we show the line between Azok and Kurdai (KURD), both sites 
on the Kazakh platform (Figure 5). The surprisingly large apparent displacement in the 
east component of this line between 1992 and 1993 was one of the most important factors 
motivating us to obtain additional measurements in 1994. Indeed, there is no resolvable 
displacement between 1993 and 1994, supporting the inference that the 1992 survey was 
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anomalous. The data displayed in this figure provide testimony to the benefit of making 
multiple occupations of sites, spaced at intervals of weeks or more early in the course of 
an experiment, in order to detect anomalies resulting form long-term systematic errors or 
blunders. And once again, the scatter of the "small" 1994 experiment is better than the 
1992 experiment, and comparable to that in 1993. 

These results were presented at the Fall 1994 DOSE meeting in Boulder. The abstract 
of the presentation is included as an appendix. 

Training in Data Analysis 

Through a special NSF initiative to help provide infrastructure support for scientific 
organizations in the former Soviet Union, PI Michael Hamburger obtained funding to 
establish two data processing facilities in Bishkek - one based at Kyrgyzgeodeziya, and a 
second facility established in collaboration with the IRIS data processing center at the 
IVTRAN laboratory. Funds were used to provide Sun Sparcstations with Gbyte disk 
drives. Exabyte tape drives, compact disk readers, laser printers and other peripherals-all 
the required hardware to support independent analysis of locally collected GPS data at 
either of these data analysis centers. 

In June, 1994 Professor Thomas Herring of MIT, assisted by MIT graduate students 
Barbara Souter and Svedana Panasyuk, conducted a formal training program on GPS data 
processing for our Russian, Kyrgyz, and Kazakh counterparts. The two-week training 
session was hosted at the IVTRAN facility, with computers provided by IVTRAN and 
Kyrgyzgeodeziya. Scientists from IVTRAN, Kyrgyzgeodeziya, Kazgeodeziya, and IPE 
Moscow participated. The training was particularly successful for the IVTRAN scien- 
tists, who were able to carry out a complete analysis of the data collected during the 1994 
field season. 


FUTURE WORK 

(1) Based on the success of last summers' field experiment, we are currently planning 
to make repeat measurements of the 86-station Tien Shan GPS network in 1995 using the 
same approach. Our counterparts in Kyrgyzstan and Kazakhstan are pursuing the neces- 
sary formalities and seeking funds to support this effort. The present plan is to send 3 
TurboRogue receivers and 6 Trimble SSE’s to Central Asia for a period of a year, where 
our colleagues will make all of the measurements. In an expansion of the observing sce- 
nario conducted last year, we will deploy the 3 TurboRogue receivers as permanent 
receivers — one at the IVTRAN polygon in the north-central Tien Shan, a second near the 
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Talgar seismic station east of Almaty (Kazakhstan), and a third near the Toktogul 
Reservoir, in the central Tien Shan. We are working on plans to make data from the 
IVTRAN receiver available via Internet. The 6 SSE's will be used to make 2 occupations 
of the network, with at least 36 hours of continuous data available from each mobile 
receiver for each site occupation. In addition, we plan to send an 2 additional SSE s for a 
few months to be used to tie in reference marks at the sites. We are exploring the use of 
an innovative new mark, developed by our colleagues at IVTRAN, which permits precise 
positioning of an antenna over a mark using a vertically positioned screw-in steel rod, and 
obviating the need for error-prone antennas and optical plummets currently in use for 
temporary field sites. 

We are also exploring possibilities of joint observations, with Chinese colleagues, that 
would permit us to tie the Tien Shan network to the Tarim Basin to the south. Because of 
the lengthy duration of our field season, we expect to be able to tie our network to other 
regional networks operating in Nepal, India, and China. These will make possible large- 
scale measurements of the India-Eurasia tectonic deformation. 

(2) As should be clear from the brief summary above, further analysis of the existing 
data and development of improved software is ongoing. We are fortunate to have 
acquired as a collaborator Tom Herring, who has used this unique data set to demonstrate 
the capability of obtaining high-precision GPS geodetic measurements even in the 
absence of local fiducial control. This very useful collaboration has developed as a result 
of the close cooperation between those who undertake the massive logistics needed to 
collect the data and those who develop the sophisticated software and expertise needed to 
process it. 

(3) At the same time, several of us are independently pursuing related projects on 
structural geology, neotectonics, seismicity, and lithospheric structure of the Tien Shan. 
These include: (a) a regional project to examine active geologic structures and seismicity 
associated with the deformed sedimentary basins in the Tien Shan (M. W. Hamburger 
and S. Ghose, Indiana University); (b) a detailed study of the surface faulting and after- 
shocks associated with the 1992 Suusamyr earthquake (R. Mellors, S. Ghose, M. W. 
Hamburger, Indiana University); (c) ongoing seismotectonic studies using the IRIS 
broadband seismic network extending across the northern Tien Shan mountain front in 
Kyrgyzstan and Kazakhstan (R. Mellors, M. W. Hamburger, G. Pavlis, Indiana 
University; through support from IRIS); (d) ongoing seismotectonic and earthquake pre- 
diction studies in the Pamir-Tien Shan region (M. W. Hamburger, G. Pavlis, Indiana 
University); (e) an interpretation of the ~ 8-mm/a, east-southeast velocity of the VLBI 
receiver near Shanghai with respect to the Eurasian plate (P. Molnar, MIT), (f) a bound 
on the Late Quaternary slip rate along the right-lateral Talaso-Ferghana fault, which is 
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based on field work carried out in 1991 with support from elsewhere (P. Molnar, MIT), 
and (g) as a synthesis of large earthquakes and geological constraints on rates of active 
deformation in the Tien Shan, a study begun with support from a previous NASA grant 
(P. Molnar, MIT). In addition, we are preparing a proposal for a multidisciplinary project 
that will include geological, geomorphic, geodetic, seismological, and magnetotelluric 
investigations in order to examine the deep structure and lithospheric dynamics of the 
Tien Shan. This project began with a workshop and field trip in the Kyrgyz Tien Shan 
during the summer of 1994. Support will be sought from the Continental Dynamics 

Program of NSF. 


TABLE 

Table 1: Site numbers, names, 4-character ID’s, the GPS days observed, and coordi- 
nates of sites in our network that were reobserved in 1994. 

FIGURES 

Figure 1. Map of central and eastern Asia showing the topographic and structural 
setting of the Tien Shan (from J.-C. Thomas, Ph.D. thesis, Univ. of Rennes, 1993). 

Figure 2. Estimates of the velocities for sites observed in our test experiment in 1994. 
The 12 sites in the northern part of the network with smaller error ellipses were observed 
in 1992, 1993, and 1994, whereas the remaining sites were observed only in 1993 and 

1994. 

Figure 3. Temporal evolution of the baseline estimates for AZOK-CHOL. Relative 
coordinates are given for 7 days in 1992, 13 days in 1993, and 9 days in 1994. The 
straight lines represent the velocity estimates from the global analysis. 

Figure 4. Temporal evolution of the baseline estimates for AZOK-MEDE. Relative 
coordinates are given for 12 days in 1992, 3 days in 1993, and 2 days in 1994. The 
straight lines represent the velocity estimates from the global analysis. 

Figure 5. Temporal evolution of the baseline estimates for AZOK-KURD. Relative 
coordinates are given for 3 days in each of 1992, 1993, and 1994. The straight lines rep- 
resent the velocity estimates from the global analysis. 
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Appendix 


Apparant Displacements From Three GPS Surveys of the Tien Shan of 
Kyrgyzstan and Kazakhstan: Tectonic Displacements or Systematic Errors? 

K Abdrakhmatov 1 , S A Aldajanov 2 , V N Gulini B H Hager 4 , M W Hamburger 6 T A 
Herring 4 , V A Kuzmichenok 2 , V I Makarov 6 , P Molnar 4 J Normandeau 7 B J 
O'Neill 7 , V Z Ostroumov 8 , S V Panasyuk 4 , M T Prilepin^, R E Reilinger 4 , I S 
Sadybakasov 9 , N A Shcherbakov 8 , B J Souter 4 , Y A Trapeznikov 11 , V E Tsurkov 12 , T 
M vanDam 4 ^Kyrgyz Inst. Seism., Bishkek, Kyrgyzstan; 2 Kazselezacshita, Almati, 
Kazakhstan; ^inst. of Physics of the Earth, Russian Acad. Sci., Moscow, Russia; 4 Dept. 
of Earth, Atmospheric, and Planetary Sciences, M.I.T., Cambridge, MA; 6 Dept of 
Geological Sciences, Indiana Univ., Bloomington, Indiana; 6 Neotectonics and Remote 
Sensing Dept., Russian Acad. Sci., Moscow, Russia; 7 UNAVCO, Boulder, CO; 
8 Kazgeodesia, Almati, Kazakhstan: 10 Kyrgyz Inst, of Geology, Bishkek, Kyrgyzstan; 
HlVTRAN, Bishkek, Kyrgyzstan; 1 2 Kyrgyzgeodesia, Bishkek, Kyrgyzstan) 

The Tien Shan (literally "Heavenly Mountains") of Central Asia comprise the 
world’s outstanding laboratory for the study of intracontinental mountain building. Under 
joint DOSE/NSF support, we are to carry out two major GPS experiments, in 1993 and in 
1995, to measure the crustal velocity field in Kyrgyzstan and Kazakhstan. In 1992, in 
preparation for our first major campaign, we established a 13-site network spanning the 
Kungey-Alatau range on the northern margin of the Tien Shan. We remeasured this net- 
woik in 1993 as pan of a larger 85-site network covering an area approximately 300 km by 
700 km and including much of the central and western Tien Shan. (This network shares 
benchmarks established by a German-Kazakh-Kyrgyz-Russian-Uzbek team in 1992.) The 
results are intriguing: Our most frequently observed baseline, extending across the 
Kungey-Alatau, from north of Issyk Kul to the Ili basin, apparently shortened by 7 ± 1 
mm, consistent with indications of horizontal shortening from geological and seismological 
studies. But there are apparent displacements of sites on the Kazakh platform, where ex- 
pected deformation rates are low, that are up to twice as large. 

In 1994, we carried out a low-budget experiment that included about half the sites 
observed in 1993. We wanted to determine if the velocities estimated from the first two 
campaigns were correct and whether a less intense experiment, with fewer receivers oper- 
ated over a longer time (with Pi’s participating in training and planning, but not in the field 
work itself), reduces cost without degrading the geodetic results. The data quality is excel- 
lent and the experiment was a great success. Our current analysis shows apparent N-S 
convergence of > 20 mm/yr on baselines spanning the mountain belt, consistent with geo- 
logic estimates. However, we find that these results are quite dependent on the data analy- 
sis scheme used (e.g. elevation cutoff, parameterization of the troposphere). In order to 
provide tectonically meaningful results, we need both a better understanding of systematic 
errors in GPS and more observations. Additional observations could be provided econom- 
ically by providing receivers for use over extended periods of time by our local 
collaborators. 



